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Abstract 

Schematic su(2) © interaction Hamiltonians, where su{2) plays 
the role of the pseudo-spin algebra of fermion operators and /13 is the 
Heisenberg algebra for bosons, are shown to be closely related to cer- 
tain nonlinear models defined on a single quantum algebra su q (2) of 
quasifermions. In particular, su q (2) analogues of the Da Providencia- 
Schiitte and extended Lipkin models are presented. The connection 
between q and the physical parameters of the fermion-boson system 
is analysed, and the integr ability properties of the interaction Hamil- 
tonians are discussed by using polynomial algebras. 



1 Introduction 

The aim of this contribution is to discuss the equivalence between systems of 
interacting fermions and bosons and systems of g-deformed effective fermions. 
The schematic fermion-boson interaction Hamiltonians that we shall deal 
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with are of two different types: 

H ++ = io f (T + u h B^B + G (T* {B^Y + T k (Bf) (1.1) 
= W/ (T + ft) + u b B*B + G (r* (5) p + T k (B*) p ) (1.2) 

and in this paper we shall restrict our study to the cases with k,m — 1,2. 
The (two-level) fermions are represented by the collective pseudospin oper- 
ators {T ,T±} and the energy difference between fermion levels is fixed by 
the energy scale Uf. Fermions are coupled to an external boson field with 
frequency Ub that is quantized through the creation (annihilation) operators 

bHb). 

Many interesting quantum models fall into this category of Hamiltonians. 
In particular, H + _ with k = p = 1 is the well-known Dicke model in Quan- 
tum Optics [1], that exhibits many of the characteristic features of quantum 
nonlinear phenomena [2] . Models based on the coupling between bi- fermions 
and bosons have been introduced long ago [3, 4] (see also the review paper [5] 
and references therein) and are particularly suitable to describe condensation 
phenomena and transitions from fermionic to bosonic phases. Among them, 
we shall study the Da Providencia-Schutte model [6] (that corresponds to the 
case H ++ with k — p — 1) and two different extensions (k = 2 and p — 1) [7] 
of the Lipkin-Meshkov-Glick model (k = 2 and p = 0) [7, 8]. 

It turns out that all the abovementioned models can be replaced by ef- 
fective su q {2) quasifermion Hamiltonians with no boson operators, provided 
the deformation parameter q is fitted in a suitable way in terms of physical 
constraints. We will refer to numerical analysis that strongly confirm this 
statement [9] and we shall comment on some of the integrability properties 
of this class of interactions. 

In this respect, we recall that the exact solvability of the su q (2) interaction 
Hamiltonian 

H q nt = q^(f + + fj)q^ (1.3) 

was already found in [10], and this operator will be used as the building block 
for the effective g-Hamiltonians that we are going to introduce. Explicitly, 
the eigenvalues of (1.3) are just the (/-numbers [2 m], (m = — j, . . . , +j), 
and the associated eigenvectors can be expressed in terms of g-Krawtchouk 
polynomials [11]. On the other hand, the use of quantum deformations of 
su(2) in this context is quite natural from a purely algebraic point of view, 
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since it is well known that fermion-boson models are related to polynomial 
generalizations of the su(2) algebra (see, for instance, [12]). In particular, if 
we consider the operators 



o — 



Tr 







K + = Tl (flt)p 



K_ = T k _ (B) p 



(1.4) 



and we take a representation of them on a certain invariant subspace of the 
interaction Hamiltonian, we get a commutation rule of the type 



where F(Kq) is a polynomial of Kq with degree (2 k + p — 1), and the same 
is true if we consider the operators K + = (B) p and K_ = (B^) p . 

In this paper, we shall use this polynomial algebra approach to get new 
(to our knowledge) integrals of the motion for all the models under study. 
Moreover, since the coefficients of the polynomial F(Kq) are given in terms of 
the physical parameters of the model (degeneracy of the fermion shells, quan- 
tum numbers of the invariant subspace, etc.) we will be able to obtain new 
interpretations of the transition to certain dynamical regimes (for instance, 
the strong field limit [2]) in terms of algebraic transformations like contrac- 
tion processes [13]. In general, we expect that the comparison between the 
algebraic properties of such polynomial algebras and those of su q (2) will ex- 
plain in more fundamental terms the efficiency of quantum algebras in order 
to model effective fermion-boson interactions. 

2 The Da Providencia-Schiitte model 

The model proposed by Da Providencia and Schutte (DPS) is a solvable 
model which exhibits a phase transition between nucleonic and pionic con- 
densates and consists of iV = 2 Q fermions moving in two single-shells each 
with degeneracy 2 Q. The DPS Hamiltonian reads [6] 



where G is the strength of the interaction and {T , T±} are the generators of 
the su(2) algebra of collective fermions: 





H = Lo f (T + u b B^B + G(T + B ] + T-B), 



(2.1) 



[T ,T + ] = T + , [T 0) r_] 



T_, [T + ,T. 



] = 2T . 



(2.2) 
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The Hamiltonian (2.1) commutes with the operator 

P = B^B- (To + Q). (2.3) 

Therefore, the matrix elements of H can be computed in a basis \mn,n) 
labeled by the eigenvalues of the number operators for fermions and bosons. 
In this basis the eigenvalues of P are given by 

P\m n ,n) — (n - m n - Q)\m n ,n). (2.4) 

In particular, we shall diagonalize H in the subspace spanned by the states 
|mn, L + rriQ + Q) = \rriQ] L, Q) which have a fixed eigenvalue L of P 

P\m n ;L,n) = L\ma;L,n). (2.5) 

In this subspace, the non-zero matrix elements of H are [9] 

(m n ; L, Q\H\m n ; L, Q) = u b L + {u f + u b )(Q + m fi ), (2.6) 

(ma + 1; L, £l\H\ m n ; L, Q) = 

Gy/(n + m n + -m n )(L + n + m n + 1). (2.7) 

The dimension of such subspace depends on the sign of L. For L > the 
quantum number ttiq can take the values 

m n = + (2.8) 

and the subspace has dimension 2Q + 1. If L < 0, rriQ takes the values 

m n = -L-n, -L-Q+ (2.9) 

and accordingly, the dimension of the invariant subspace is 2 Q + L + 1. 

In general, we shall say that the system is in a normal phase when the 
correlated ground state is the eigenstate of the symmetry operator P with 
the eigenvalue L = 0. The denomination deformed phase will be assigned to 
cases where the correlated ground state is an eigenstate of P with eigenvalue 
L 0: if L > we shall have a so-called bosonic phase and the case L < 
corresponds to a fermionic one. 
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2.1 The DPS algebra 

Let us define the operators 

K = T K + = T + B ] K_ = T_ B. (2.10) 

For a given Q and L > 0, if we consider the action of these operators within 
the subspace |toq; L,Q) the following DPS algebra is obtained: 

[K ,K ± ]=±K ± 

[K + ,KJ [ = -n(n+l) + {l + 2L + 2Q)K + 3Kl (2.11) 

It can also be checked that the DPS algebra for the L < case is just (2.11) 
where L is replaced by —L. Therefore, a new integral of the motion for the 
DPS model is given by the Casimir operator of the DPS-algebra, that can 
be found by standard methods [14] and reads 

C = K + + Kl + ( L + n - 1) Kl - {( L + Q) + Q(Q + 1)} K . (2.12) 

In the corresponding L > subspace, the eigenvalue of C is 

n(n + i)(n + L). (2.13) 

We remark that the DPS algebra (2.11) is a quadratic generalization 
of the su(2) Lie algebra. In fact, the su(2) algebra can be obtained as a 
contraction [13] of the DPS algebra in the L — > oo limit. Namely, if we 
define the "contracted generators" as 

J = K J ± = -j=K± (2.14) 

and we compute their commutation rules we get 
[J ,J±]=±J ± 

[J+,J-] = y{-fi(fi+ 1) + (1 + 2L + 2fi) J + 3J 2 } (2.15) 

which is still isomorphic to the DPS algebra. However, the L — > oo limit of 
these commutation rules gives 

[J ,J ± ]=±J ± [J + ,J_] = 2J . (2.16) 
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In physical terms, the L — > oo limit is just the well-known "strong field" 
limit of the Dicke model in Quantum Optics [2], for which the interaction 
dynamics is given by su(2). If we rewrite the DPS Casimir in terms of the 
contracted generators we find 

C = L J + J_ + J 3 + ( L + ft - 1) 4 - {( L + ft) + ft(ft + 1)} J (2.17) 

and the su(2) Casimir operator C su ( 2 ) is obtained by computing 

C su{2) = lim % = J+ J_ + J (Jo - 1). (2.18) 

L-^oo J_j 

Note that the contracted eigenvalues are just ft(ft + 1), as it should be. As 
we shall see in the sequel, this polynomial algebra approach can be applied 
to all the fermion-boson Hamiltonians under consideration. 



2.2 Effective su q (2) Hamiltonians for the DPS model 

The quantum algebra su q {2) is a Hopf algebra deformation of su(2) [15] with 
generators {T ± ,T } and commutation rules 

[f ,f ± ] = ±f±, [f + ,f_] = [2T ]„ (2.19) 

where the g-number [x] q is defined by 

[x\ q = qX ~ g "' = Sinh ^f (2.20) 
q — q~ L sinh(z) 

We shall use alternatively q and z (where q = e z ) as the deformation param- 
eter, and we shall assume that q is real. The su{2) algebra (2.2) is recovered 
from (2.19) in the limit q -> 1 (z -> 0). 

When q is not a root of unity, the irreducible representations of su q {2) 
are obtained as a straightforward generalization of those of su{2) [16, 17]: 

T \j,m) = m \j,m), 

T+ \j, m) = yf\j + m + l] q \j -m] q \j, m + 1), 

T-\j, m) = - m + l] q [j + m} q \j, m-l). (2.21) 
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By following [9, 10], we consider an effective Hamiltonian defined as 

H q = co b L + (uj b + uf)(f + + X (q)q^(T + + f_)q^ (2.22) 

where x(o) * s a scalar function to be fixed, and H q will be realized in a 
su q (2) irreducible representation with the same dimension as the subspace 
spanned by |toq;L, Q) (therefore, with j = j(Q,L)). The non-vanishing 
matrix elements of (2.22) read 

(j, m\H q \j, m) = u b L + (ujf + uo b )(m + Q), (2.23) 
(j,m+l\H q \j,m) = x(q)q {m+ ^^/[j+m + l] q [j-m} q . (2.24) 

In order to fit the dimension of the representation with respect to the invari- 
ant subspace of the DPS model, we have to take j — VI and m = rrin for the 
effective L > model, while for L < 0, j = Vl + |f and m = run + \- 

The main conclusion of [9] (see also [10] for the Dicke model) is that 
the Hamiltonian (2.1) is essentially equivalent to (2.22). In other words, the 
bosonic degrees of freedom included in (2.1) may be reabsorbed by the q- 
deformation in (2.22) provided that q is defined as an appropriate function 
of both Q and L. In this way it is possible to regard H q as an effective Hamil- 
tonian with physical properties similar to those of H. As it is extensively 
shown in [9] through numerical studies, both the ground state energy and 
the full spectrum of the DPS model can accurately be reproduced by using 
the effective (/-Hamiltonian (2.22). 

2.3 The q-DPS algebra 

By following the same algebraic approach leading to the previous DPS alge- 
bra, now we should consider the su q {2) operators 

K = f K + = q^ T+ q^ K_ = q^ T_ q^ (2.25) 

such that the effective Hamiltonian (2.22) is a linear function of K and K±. 
In this new basis, the ^-commutation rules of su q {2) read 

[K , K ± ] =±K ± , q K + X_ - q- 1 K_ K + = q 2K ° [2 K ] q (2.26) 
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and these expressions hold for any irreducible representation j of su q (2). The 
Casimir element for this algebra is 

C q = [K ] q [K - l] q + q- 2Ko+1 K + K_, (2.27) 

and its eigenvalue is just [j] q [j + l] q . Obviously, C q is an integral of the 
motion for H q (2.22). By working on a fixed irreducible representation j, 
the latter g-commutator can be rewritten as the following commutation rule 
(that hereafter we will call the g-DPS algebra): 

[K + ,K_] = q +1 ^ -\^- 2 q iK \ (2.28) 
1 — q z 1 — q z 

which should have algebraic properties closely related to the ones of the 
DPS algebra (2.11), since both models are physically equivalent. In fact, 
the analytic fitting q = g(ft, L) should be found by comparing the properties 
of the DPS and g-DPS algebras. Work on this open problem is actually in 
progress [22]. 



3 Extended Lipkin models 

As a second example of Hamiltonians including fermionic and bosonic degrees 
of freedom, let us introduce the extended Lipkin-Meshkov-Glick Hamiltonian 
(LE model) [7], which is just (1.2) with k = 2 and p — 1: 

H+_ = u f (T + ft) + u b B*B + G{TlB + T 2 B ] ). (3.1) 

which commutes with the operator [18] 

P (+) = BiB + ±(T + n). (3.2) 

Therefore, the matrix elements of the LE model can be calculated in a 
basis labeled by the eigenvalues of P(+y 

P(+)|ran,n) = L\m n ,n) = (n + ^(ft + m n ))\m n , n), (3.3) 
and we shall consider the invariant subspace with L fixed: 

\msi,L - i(ft + ma)) = \m n ; L, ft). (3.4) 
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In this subspace the non-zero matrix elements of H read 

(mn] L, Q\H + -\mn; L, Q) = cu b L + (cuf - ^u b )(Q + mn), (3.5) 
(m n + 2; L,tt\H + -\m n ; L,Q) = G\ L - htt + mn) x 



'(n + m n + 2)(n + m n + l)(n - m n )(0 -m n - 1), (3.6) 

and we have to distinguish the following classes of invariant subspaces: 

L> £l,L integer, m n + = 0, 2, 20, 

L > Q, L half integer, mo + = 1,3,..., 20 — 1, 

L < 0, L integer, m n + = 0, 2, 2L, (3.7) 

L < O, L half integer, mn + O = 0, 2, 2L — 1. 

Another Lipkin-type Hamiltonian can also be defined through (1.1) with 
k = 2 and p = 1: 

#++ = W/ (T + O) + a^fl + G(T%Bt + T 2 ^), (3.8) 

which differs from (3.1) in the ground state correlations [7]. Since H ++ 
commutes with the operator 



(-) 



B^B--(T + n), (3.9) 



its matrix elements can be calculated in a basis labeled by the eigenvalues of 
P(_). Namely 

P(_)|m n ,n) = (n - ^(O + m Q ))|mn,n), (3.10) 

and in this case we shall compute the matrix elements within the subspace 
spanned by the states \mn,L + |(0 + mn)) = |m^;L, O). The non-zero 
matrix elements of H ++ (3.8) are now given by 

(mn] L, Q\H ++ \m n ; L, O) = + (a;/ + ^u; 6 )(fi + mn), (3.11) 
(m n + 2; L, Q\H ++ \m n ; L, O) = CWL + J(Q + m n ) + 1 x 



(Q + m n + 2){Q + mn + - m n )(0 - m n - 1), (3.12) 
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where the dimension of the subspace depends on L and f2, since we have 
to consider the following possibilities for the set of values of the quantum 
number ttiq: 



L > 0, L integer, + Q — 0, 2, 2f2, 

L > 0, L half integer, m n + Q = 1, 3, 2Q - 1, (3.13) 

L < 0, L integer, m Q + = -2L, -2L + 2, 2fi, 

L < 0, L half integer, m n + tt = -2L, -2L + 2, 2Q - 1. 



3.1 LE algebras and their Casimir operators 

The polynomial algebra approach can also be used for the LE models (3.1) 
and (3.8). We start our analysis by recalling the cubic algebra linked to the 
original Lipkin-Meshkov-Glick (LMG) Hamiltonian [8], since the latter will 
appear as the strong field limit for the extended LE models. 



3.1.1 The LMG algebra 

We recall the su(2) LMG Hamiltonian [8] given by 

H = uj f T + x{T 2 + + T 2 _). (3.14) 
If we define (see also [19]) 

K = T K + = T% K_ = T 2 (3.15) 

we get the cubic algebra 

[K ,K±\ = ±2K± 

[K+, K_] = 4{2 + 1) - 1} K - 8 Kl (3.16) 

where we have identified j = Q > 0. Note that this algebra is isomorphic to 
the Higgs algebra [20, 21] for any value of Q. The Casimir operator for this 
algebra is found to be: 

C = K + — Kq+4Kq + {2n(Vl + 1) - 5} Kq - 2{2Vl(n + 1)-1}K (3.17) 

and the eigenvalues of this operator are f2(f2 2 — l)(fi + 2). 
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3.1.2 The LE+ algebra 

From (3.1), we can consider the operators 

K = T K + = TlB K_ = T 2 B ] . (3.18) 

If we compute its action on an invariant subspace of the type L > Q, we 
obtain the following quartic generalization of the su(2) Lie algebra: 

[K ,K±] =±2K ± 

[K+, K_\ = a + a 1 K + a 2 K 2 + a 3 + a 4 (3.19) 

whose structure constants aii(Q,L) read 

a = tt(tt 2 - l)(n + 2) 

ai = 2(1 - 2fi(fi + -2L-1) 

a 2 = 7-6fi(fi + l) (3.20) 

a 3 = 4(ft-2L-l) 
«4 = 5. 

The Casimir operator for (3.18) can also be computed and gives a new integral 
of the motion for H + : 

C = K + K„ + /3 + faK + foKl + (3 3 K 3 + faK* + /3 5 K 5 (3.21) 
with = -tt(tt 2 - l)(n + 2) 

1 = 2(L + 1) + i^(-12 + tt(tt 2 + 6fi - 5) - 8(fi + 1)L) 
/3 2 = -6 - 5L + ifi(13 + 2fi(3 - ft) + A(n + 1)L) 

= % - n(n + 3) + 4L 

/? 4 = iQ - 3 - L 

The eigenvalue of C is found to be ~\Vt(Vt 2 - + 2){Q - 2L). It can 
also be checked that the transformation (2.14) leads to the LMG algebra 
(3.16) and Casimir (3.17) as the "strong field" contraction L — > oo of (3.19) 
and (3.21), respectively. Similar quartic algebras can be obtained for the 
remaining invariant subspaces 
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3.1.3 The LE++ algebra 

For the second class of LE models (3.8) we define the generators 

K = T K + = rj B ] K„ = T 2 _ B. (3.22) 

The associated quartic algebra for L > reads: 
[K ,K ± ] = ±2K ± 

[K + , K_\ = a + a 1 K + a 2 K 2 + a 3 + a A K* (3.23) 

with a = -n(Q 2 - l)(n + 2) 

ai = -2(1 - 2tt(tt + + 2L + 1) 

a 2 = -7 + Qtt(tt + 1) 
a 3 = -4(Q + 2L + 1) 
«4 = —5 

The Casimir operator is now 

C = K + K_ + (3 + f3 l K + f3 2 K 2 + p 3 K 3 +f3 4 K* + (3 5 K 5 (3.24) 

with fa = tt(tt 2 - l)(n + 2) 

(3 1 = 2L- ±fi(-4 + fi(fi 2 + 6fi + 3) + 8(fi + 1)L) 
/3 2 = 1 - 5L + iQ(-9 + 2il(tt - 1) + 4(Q + 1)L) 
/3 3 = -§ + ft(ft + 3)+4L 
/3 4 = -iQ + 2-L 

and the eigenvalue of C is \Q{yi 2 — l)(f2 + 2)(2 + f2 + 2L). Once more, the 
"strong field" contraction of this algebra gives rise to the same LMG algebra 
(3.16), that underlies the asymptotic L — > oo dynamics of both LE models. 

3.2 Effective su q (2) Hamiltonians for the LE models 

Like in the case of the DPS model, an effective g-Hamiltonian for (3.1) has 
been introduced in [9] 

H q = u b L + (cu f - l -uj b ){f Q + + X (q)q fo (fl + f 2 )q f °, (3.25) 
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which has the following non-vanishing matrix elements 

(j, m\H q \j, m) = u b L + (u f - ^w 6 )(Q + m), (3.26) 

(j, m + 2\ H\j,m} = x(q)q 2{m+1) x 

y/\j + m + 2] q [j + m+ l] q [j - m] q [j - m - l] q . (3.27) 

As the first step in order to fit the dimensions of (3.1) and of (3.25) the 
appropriate relation j = j(Q,L) has to be found and, as a consequence, 
m = m(mn, f2, L) [9]. Afterwards, we have to consider as the effective Hamil- 
tonian the restriction of the matrix elements (3.26) and (3.27) to the invariant 
subspace spanned by \j, m) with m = —j, —j + 2,...,j — 2,j. In this way 
we obtain the effective matrix elements 

(j,m + 2\H q \j,m) = X {q)h{L,Q,m n ). (3.28) 

For values of L > Q (L integer), we find j — Q, m — rrin and the function 
h(L, Q,mn) is (the function x(q) is also explicitly defined in [9]): 

h(L,n,rrm) = g 2(mn+1 V[fi + m n + 2] q [Q + m n + l] q [Q - m n ] q [Q - m n - l] q . 

(3.29) 

Afterwards, we have followed the same procedure as in the g-DPS model, 
and we have looked for values of the q parameter (and, consequently, of 
x(q)) which may absorb bosonic degrees of freedom of (3.1) and give rise to 
a similar spectrum for the purely fermionic g-deformed Hamiltonian (3.25). 
The second type of LE model (3.8) can also be approximated by the same type 
of g-Hamiltonian (3.25) and, in both cases, numerical computations lead to 
an excellent fitting between the LE models and the effective g-Hamiltonians 
[9]. Therefore, we conclude that certain interactions between fermions and 
bosons can accurately be described by using g-fermions as quasiparticles (i.e; 
effective fermionic degrees of freedom) under the exactly solvable interaction 
given by the Hamiltonian H q n ^ (1.3). 

From the mathematical point of view, a very interesting question to be 
solved is to find a suitable analytical expression of the deformation parameter 
q in terms of the representation space labels Q and L. In this respect, the 
quartic LE algebras previously introduced should be relevant, since the su q {2) 
model (3.25) leads to the following natural definition of the K operators 

K = f K + = q fo f\ q fo = q f ° fl q f ° (3.30) 
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and their commutation rule [K + ,K_] (which is a generalization of (2.28)) 
has to carry essentially the same algebraic information as (3.19) and (3.23). 
Another intteresting feature appears in the analysis of the effective g-Hamil- 
tonian for H ++ (3.8), since a su(2) symmetry of the model can dynamically 
be restored for certain negative values of L [9]. A complete study of the 
LE algebras together with the abovementioned algebraic properties of these 
fermion-boson interactions will be addressed in a forthcoming paper [22]. 
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